Temperature dependent band structure of the Kondo insulator 
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We present a Qantum Monte Carlo (QMC) study of the temperature dependent dynamics of 
the Kondo insulator. Working at the so-called symmetrical point allows to perform minus-sign free 
QMC simulations and thus reach temperatures of less than 1% of the conduction electron bandwidth. 
Study of the temperature dependence of the single particle Green's function and dynamical spin 
correlation function shows a surprisingly intricate low temperature band structure and gives evidence 
for two characteristic temperatures, which we identify with the Kondo and coherence temperature, 
respectively. In particular, the data show a temperature induced metal-insulator transition at the 
coherence temperature. 
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The theoretical description of the Kondo lattice re- 
mains an outstanding problem of solid state physics. This 
model, or variations of it, may be viewed as the appropri- 
ate one for understanding such intensively investigated 
classes of materials as the heavy electron metals 
and the Kondo insulators || . Experimental results indi- 
cate that the electronic structures of Kondo lattice com- 
pounds undergo quite dramatic changes with tempera- 
ture 0]. It is the purpose of the present manuscript to 
report a QMC study of the electronic structure of the so- 
called Kondo insulator, which shows that this model in- 
deed undergoes a quite profound change of its unexpect- 
edly intricate band structure as temperature increases. 
We are using a one dimensional (ID) 'tight-binding ver- 
sion' of the model with L unit cells and 2 orbitals/unit 
cell: 
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Here cj a (fj a ) creates a conduction electron (/-electron) 

in cell i, ni t(T =fJ f t a . Throughout we consider the case 
of 'half- filling' i.e. two electrons /unit cell and, as an 
important technical point, we restrict ourselves to the 
symmetric case, e/=f//2. The latter choice, while prob- 
ably not leading to any qualitative change as compared 
to other ratios of ef/U, has the crucial advantage that 
at half-filling the model acquires particle-hole symme- 
try, i.e. the Hamiltonian becomes invariant under the 
transformation aj, CT — ► exp(iQ ■ Ri)a\ a , where a—c,f 
and Q=(ir, ir, ...) (this holds for bipartite lattices with 
only nearest neighbor hopping). Particle- hole symme- 
try in turn implies that the QMC-procedure does not 
suffer from the notorious 'minus-sign problem' anymore, 
so that reliable simulations for temperatures as low as 
/3t=30, corresponding to w 0.8% of the conduction elec- 
tron bandwidth, can be performed without problems. 
This allows to scan the dynamical correlation functions 



of the system over a wide temperature range. Previouslt 
a QMC study for the 2 dimensional model at half-filling 
was performed by Vekic at al. || , more recently a study 
of the temperature dependence of static susceptibilities 
for the strong coupling version of the model has been re- 
ported by Shibata et al. 0. 

It is widely believed that the Kondo lattice has two dis- 
tinct characteristic temperatures. At the Kondo temper- 
ature, Tk, the /-electrons start to form loosely bound 
singlets with the conduction electrons. This manifests 
itself in a deviation of the spin susceptibility from the 
high-temperature Curie form due to the 'quenching' of 
the /-electron magnetic moment and an increase of the 
dc-resistivity due to resonant scattering from the newly 
formed low energy bound states Qj. The second (and 
lower) characteristic temperature is the coherence tem- 
perature, T co h, where the local singlets establish long 
range coherence amongst themselves so as to participate 
in the quasiparticle bands of a Fermi-liquid like electronic 
state. Experimentally the coherence temperature is sig- 
naled by the onset of a decrease of the dc resistivity with 
temperature |Q, and the formation of the 'heavy bands' 
which (judging by the volume of the Fermi surface) in- 
corporate the /-electrons . While there is as yet no ex- 
perimental proof, one might expect on the basis of these 
considerations, that at temperatures above T co h the /- 
electrons do not participate in the Fermi surface volume, 
whereas they do so below. 

Turning to the Kondo insulator we note that the electron 
count for these systems is such that the 'Fermi surface' 
comprising both, conduction and /-electrons, would pre- 
cisely fill the Brillouin zone so that the system is a 'nomi- 
nal' band insulator. If increasing temperature causes the 
/-electrons to 'drop out' of the Fermi surface volume, this 
should manifest itself as an insulator-to-metal transition 
because the volume of the collapsed Fermi surface is no 
longer sufficient to cover the entire Brillouin zone. Tran- 
ferring the above scenario for heavy Fermion metals to 
the Kondo insulator one would therefore expect that the 
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system remains a metal above T co h, with a Fermi surface 
that excludes the /-electrons, and becomes insulating be- 
low T co h, when the /-electrons participate in the Fermi 
surface volume to turn the system into a 'nominal' band 
insulator. In fact insulator-to- metal transitions which 
are induced by temperature 0, or hydrostatic pressure 
H have been observed experimentally. As will be seen 
below, our data for the Kondo insulator are remarkably 
consistent with such an interpretation. 
To begin with, we consider the T-dependence of the sin- 
gle particle spectral function. This is defined as (a=c, /) 



A a (k,oj) = -^e^Ki/la^l/x)! 2 ^-^), 



(2) 



where the sum is over all eigenstates \v) of H — fiN in the 
grand canonical ensemble, Z denotes the partition func- 
tion and u> ufl the difference of the energies of the states 
v and /i. Figure [l] shows the angle-integrated spectral 
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FIG. 1. Angle integrated spectral density D(oj) for an 
L—16 Kondo lattice for different temperatures. Parameter 
values are U/t=8, V/t=l. 

density, D a (uj) = (l/L)^2 k A a (k,u}). At the lowest tem- 
perature, /3t=30, D a (uj) shows the behaviour expected 
for a Kondo insulator: the c-electron density is roughly 
consistent with the standard ID tight-binding density of 
states, with indications of the two van-Hove singularities 
at ±2t. Around however, the spectral density shows a 
small but unambiguous gap, which demonstrates the in- 
sulating nature of the ground state. The /-like spectral 
density shows very sharp low energy peaks at the edges of 
this gap, as well as high-intensity 'Hubbard bands' at ap- 
proximately ±Z7 /2. There are also two weak 'side bands' 
at approximately ±0.8t. With the exception of these, all 
features are consistent with exact diagonalization 0] and 
d — > oo results jlO| at T=0. This suggests that the side- 
bands are already an effect of the finite temperature and 
the further development with T confirms this. Inspec- 
tion of the series (3t=30, 10 shows that with increasing T 



a transfer of spectral weight from the low energy peaks 
at the gap-edges into the sidebands is taking place. This 
is accompanied by a narrowing of the gap, and at /3i=10 
the gap practically closes; in the spectrum for this tem- 
perature the two low-energy peaks seem to have collapsed 
into a single one right at fj, - we believe that these are in 
fact still two peaks, which however are too closely spaced 
to be resolved. Increasing T even further (/3t—5), the 
low energy peaks disappear completely. The c-electron 
spectral density no longer shows any indication of a gap, 
whereas the /-like sidebands stay at a relatively high en- 
ergy away from /i. The extreme low energy states thus 
have pure c-character (within the resolution of the QMC 
procedure) and the system is a metal with a c-like Fermi 
surface. We therefore interpret the temperature where 
the gap closes as the analogue of the coherence temper- 
ature. With increasing T the energy of the sidebands is 
lowered and at f3t=3 they collapse into a single peak right 
at ji. This probably indicates a second transition and 
at the relatively large value of (3t=l another reconstruc- 
tion of the band structure has taken place, namely the 
disappearance of the /-like sidebands. The upper and 
lower Hubbard band for the /-electrons are now quite 
broad, and actually the possibility that there are very 
low-intensity /-like features near [i cannot be completely 
ruled out. However, the overall trend is quite obviously 
a strong decrease of the sidebands. Interpreting the lat- 
ter as the lattice-analogue of the Kondo-resonance in the 
impurity case (TTJ , the temperature of the second transi- 
tion should correspond to the Kondo temperature Tk ■ In 
the present case Tk is very high because of the relatively 
large value of the c-f hybridization, V=t. We also note 
that the closing of a gap in the /-like density of states in 
2 dimensions was previously found by Vekic et al. j5j. 
To get a more detailed picture, we consider the fc-resolved 
single particle spectral function, shown in Figure ^| for 
some temperatures in between the two transitions. For 
flt=30 the c-electrons show a standard cos(fc) band, al- 
beit with a clear gap at k F =ir/2, the Fermi momentum 
of nonhybridizcd conduction electrons. At this momen- 
tum the band changes its spectral character and 'bends 
over' into a practically dispersionless /-like band, which 
can be followed up to k—n. This kind of band structure 
is familiar from various studies §,0,H. The weak and 
practically dispersionless /-like sidebands are at some- 
what higher energy and, at very high energies, the /-like 
Hubbard bands. The width of the Hubbard bands seems 
to depend strongly on momentum - this is a deficiency of 
the QMC and maximum entropy method, which is most 
accurate near \i. As seen in the fc-integrated spectra, 
raising the temperature leads to a transfer of spectral 
weight from the 'flat band' forming the single-particle 
gap into the Kondo resonance-like sidebands. At /3i=5, 
where the gap has closed, the c-electron spectrum shows 
a very conventional cos(fc)-band with no more indication 
of any gap. At k F there is now one symmetric and un- 
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split peak right at \x - as it is required by particle-hole 
symmetry for a metallic system. The system thus has a 
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FIG. 2. Momentum resolved single-particle spectral func- 
tion A(k, <jj) at different temperatures. Parameters as in Fig- 
ure jj], the momentum k increases in steps of 2n / L—tt /8 from 
the bottom of each panel. 

Fermi surface as expected for unhybridized c-electrons, 
i.e. the /-electrons indeed have 'dropped out' of the 
Fermi surface volume. The dispersionless /-like side- 
bands are at a relatively high energy. The tiny low en- 
ergy 'foot' seen in some of the c-like spectra may indicate 
a very weak mixing of the c-electrons into the Kondo- 
resonance but apparently this no longer leads to a gap. 
Rather, the Kondo resonance is now essentially decou- 
pled from the Fermi surface physics. Finally, for /3t—l, 
even these Kondo resonance-like sidebands have disap- 
peared and the only /-like peaks in the spectral func- 
tion are the upper and lower Hubbard bands. At this 
high temperature the /-electrons do not participate in 
the low energy physics at all. The expectation value of 
~VJ2io-( c l.crfia + H.c.) decreases by as 30% between 
/3i=30 and (3t—l - while there is appreciable mixing even 
at high temperature, this does obviously not lead to co- 
herent band formation any more. 

We proceed to the dynamical spin correlation function 
(SCF), defined as 



S(k,u) = i^e-^|H^(fe)| M )| 2 ^-^), 

U,fJ, 



(3) 



component of the spin-operator for a-electrons. This is 
shown in Figure [| At /?i=30, the /-like SCF shows an 
intense branch of low energy excitations with a tiny but 
clearly resolved spin wave-like dispersion. The spectral 
weight of this branch is sharply peaked at k—ir, indi- 
cating relatively long ranged and strong antiferromag- 
netic spin correlations. Fitting the equal (imaginary) 
time /-like spin correlation function in real space to the 
expression S(r)=A(e~ r ^+e~( L ~ r ')/ ( ') we obtain the val- 
ues (=4.61,4.67 and 2.18 for /3t=30, 20 and 10. The 
dominant feature in the c-like SCF on the other hand, 
is a free electron-like particle-hole continuum. Interest- 
ingly enough, there is also a replica of the /-electron spin 
wave branch in the c-like SCF. This shows that at low 
excitation energies c and / electrons behave as a single 
'all-electron fluid', whose excitations have composite f-c 
character. The free electron continuum itself does have a 
gap of w 0.6i at k=ir - this corresponds to approximately 
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FIG. 3. Dynamical spin correlation function S(k,uj) for dif- 
ferent temperatures. Parameters are as in Figure hi 

twice the single particle-gap in A{k,ui) We pro- 

ceed to f3t=5, where the single-particle gap has closed. 
The /-like SCF still shows a low-energy peak with fi- 
nite excitation energy, which however is practically dis- 
persionless, both with respect to its energy and with re- 
spect to its spectral weight. In other words, the magnetic 
/-excitation becomes practically immobile. The metal- 
insulator transition thus is related to a drastic change of 
the spin correlation function p| - the question whether 
the longer ranged spin correlations below the transition 
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are the driving force behind the gap formation |13(] or 
whether the change of the spin correlations is merely 
a 'byproduct' of the collapse of the single-particle gap, 
remains to be clarified. In any way, the almost com- 
pletely localized spin dynamics of the /-electrons natu- 
rally should lead to a Curie-law for the static spin suscep- 
tibility at temperatures above T co h- A crossover from an 
activation-gap dominated susceptibility at low tempera- 
tures and a Curie law at high temperatures has indeed 
been observed by Shibata et al. In the c-like SCF 
the particle-hole continuum persists and the gap near 
k—TT is now very small or zero (the absence of the gap 
in the single particle spectrum suggests it to be zero). 
There is no more distinguishable peak in the c-like SCF 
which would correspond to the dispersionless /-like spin 
resonance - this suggests that the /-electrons now are 
largely decoupled from the c-like band, as indicated by 
their non-participation in the Fermi surface. At the very 
high [3t=l, there is still some (very weak) indication of 
the low energy /-electron spin resonance, but the inten- 
sity is low and the resonance is now relatively broad. It 
should be noted that the relative change of the /-like 
magnetic moment is less than 2% over the entire temper- 
ature range we studied - temperature thus affects only 
the coupling of these moments to the conduction elec- 
trons. 

For the relatively large value of V=l the higher of the two 
crossover temperatures (i.e. the Kondo temperature) is 
already rather high, (3t>l. Based on the impurity results 
JLlj one might expect that for smaller V the Kondo tem- 
perature is lower, and to check this, we have performed 



spin correlation function shows no more indication of the 
low energy resonance - the overall picture is completely 
the same as for V/t=\ and /3t=l, with the sole excep- 
tion that all features are much sharper due to the lower 
temperature. Here we do not pursue the issue of the de- 
tailed parameter dependence of the characteristsic tem- 
peratures - it is quite obvious, however, that lower values 
of V shift the characteristic temperatures of the system 
towards lower values, but otherwise leave the physics un- 
changed. 

In summary, we have studied the temperature evolution 
of various dynamical correlation functions of the Kondo 
insulator and found indications for two distinct electronic 
crossovers. At the low temperature crossover the single- 
particle gap closes so that the system becomes a metal, 
the magnetic correlations on the /-sites become local- 
ized, and the 'spin gap' closes due to c-like spin excita- 
tions. While c and / electrons seem to form a coherent 
'all-electron fluid' below the crossover temperature, the c 
and /-like features in the correlation functions above this 
temperature are decoupled. We therefore interpret this 
temperature as the analogue of the coherence tempera- 
ture in heavy-Fermion metals. At the high-temperature 
crossover both, the dispersionless /-like Kondo-resonance 
in the single-particle spectrum and the /-like low-energy 
spin excitation disappear. The only remaining /-like fea- 
ture in the single particle spectrum are the high-energy 
'Hubbard bands', corresponding to the 'undressed' tran- 
sitions f 1 — > /° and f 1 — > f 2 . We therefore interpret 
this second temperature as the Kondo temperature of 
the system. 
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FIG. 4. Single particle spectral density D(ui) at /3t=5 for 
different V/t. All other parameters are as in Figure |l]. 

simulations at fixed /3t=5, but with variable hybridiza- 
tion. Figure |] shows the results for the single particle 
spectral density. The weight of the Kondo resonance- 
like sidebands for fixed temperature decreases with V, 
and for V/t=0.25 they have disappeared completely. In 
the A;-resolved c-electron spectrum (not shown) there is 
now only a very sharp nearest neighbor-hopping band 
with a clear Fermi level crossing at k° F and the /-like 
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